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METHOD FOR CONTINUOUS SWEEPING AND 
SEPARATION OF MULTIPLE SEISMIC VIBRATORS 



[0001] This application claims tiie benefit of U.S. Provisional Patent 
5 Application No. 60/494.194 filed on August 1 1 . 2003. 

FIELD OF THE INVENTION 

[0002] This invention relates to the field of seismic data acquisition. 
Specifically, this invention is a method of acquiring seismic data using multiple 
vibrators without a multiple sweep listening time requirement. 

10 BACKGROUND OF THE INVENTION 

[0003] Seismic vibrators have long been used in Vne seismic data 
acquisition Industry to generate the acoustic signals needed in geophysical 
exploration. The conventional use of vibrators involves several well- 
understood steps. First, one or more vibrators are located at a source point 

15 on the surface of the earth. Second, the vibrators are activated for several 
seconds, typically ranging from four to sixteen, with a pilot signal. The pilot 
signal is typically a sweep signal that varies in frequency during the period of 
time in which the vibrators are activated. Third, seismic receivers are used to 
receive and record response data for a period of time equal to the sweep time 

20 plus a listen time. The period of time over which data is recorded Includes at 
a minimum the time necessary for the seismic signals to travel to and reflect 
off of the target reflectors of interest, and for the reflected signals to return to 
the receivers. Fourth, seismograms are generated by cross con-elating the 
recorded data with either the pilot signal or a reference sweep. Fifth, the 

25 sweep and correlation steps are repeated several times, typically four to eight, 
and the conrelations are added together in a process referred to as stacldng. 
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Finally, the vibrators are moved to a new source point and the entire process 
is repeated. 

[0004] Several problems are known to exist with conventional vibrator 
technology. First, the con-elation process is Icnown to result in oonrelatlon side 
5 lobes, which can Influence the accuracy of the final processed data. Second, 
vibrator hamnonic distortion results in noise, known as hamnonlc ghosts, after 
conreiation with the pilot A partial solution to this problem is the use of 
upsweeping pilot signals, in which the sweep starts at low frequencies and 
increases to high frequencies. This approach places the correlation ghosts 

10 before the main correlation peak where they will not Interfere with later, and 
hence weaker, reflections. In addition, to minimize noise from hanmonics, 
multiple sweeps are perfomied with incremental phase rotation of the sweeps 
so that after correlation and stack, the harmonics are reduced. For example, 
to suppress harmonics through fourth order, four sweeps might be perfomned 

15 with a phase rotation of 360 degrees divided by four, i.e., 0, 90, 180, and 270 
degrees. The data are stacked after correlation with harmonics accordingly 
reduced, although not eliminated. Third, in order to accurately process the 
recorded data, both the sweep time and a listen time must be included in the 
recording time of the seismic receivere for each sweep. The listen time is 

20 Important to ensure that the resulting data from each sweep can be accurately 
processed. In addition, multiple sweeps are often required to inject sufficient 
energy Into the ground. Multiple short sweeps can result in better data quality 
than long sweeps through the use of phase rotations to reduce hamnonic 
noise and by reducing ground roll reverberations. However, the use of 

25 multiple sweeps with each sweep followed by a listening time limits the rate at 
which energy can be put into the ground and the survey acquired. Fourth, the 
recording of high frequencies can be limited by the simultaneous recording of 
the signals from an array of vibratore, each vibrator at a different position and 
elevation and having a different coupling with the ground. 

30 [0005] The cost of land sun/eys largely depends on the time it takes to 
record the sun/ey, and cost is affected by the length of time required to record 
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data at each source station as well as the time It takes to move the vibrators 
to the next station. The time taken to record data at each source station 
depends on the number of sweeps, the sweep length, and the listen time. For 
example. If four 8-second sweeps are performed, each having a 7-second 

S listening time, at least 60 seconds Is required at each station. Typical data 
acquisition systems also require 3-5 seconds before they are ready to start a 
new record, which can add another 12-20 seconds to the time at the source 
station. If multiple stations could be recorded simultaneously, or the need for 
a listening time reduced or eliminated, then less time would be needed for 

10 recording the survey, therefore reducing the overall cost of the survey. 
Similarly, improved methods of reducing side lobe correlations and harmonic 
ghosts would improve the quality of land surveys. 

[0006] In 1995. Andersen in U.S. Patent No. 5,410,517 disclosed a 
method to cascade vibrator sweeps to eliminate unproductive listening times. 

IS while still maintaining the advantages of using multiple short sweeps. The 
metiiod includes an incremental phase rotation of the sweep segments and 
the use of a second cascaded sweep with an extra sweep segment to 
suppress hamnonic ghosts after correlation with the reference. For example, 
but without limitation, in a seismic acquisition program in which it was desired 

20 to reduce up to tiie fourth order hannonic, four sweep segments with an 
appropriate phase rotation are required. The phase rotation angles could be 
0, 90. 180, and 270 degrees, respectively, although other choices could also 
be made. If eight-second sweep segments were used and 7 seconds 
listening time is required, then the total record time is 39 seconds. This 

25 compares to a standard sweep and listen time for 4 sweeps, which is 60 
seconds. Combined with Andersen's prior U.S. Patent No. 5,347,494, which 
disclosed a method of producing simple seismic wavelet shapes with minimal 
side lobe energy, improved quality vibrator data can be obtained. 
Nevertheless, even with this improved quality data the limitations of the 

30 con'elation pn3cess, problems with hamnonlcs ghosting, and array limitations 
are present. 
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[0007] Another method used by Industry to increase the rate of seismic 
acquisition is to use more than one vibrator and record multiple source 
locations simultaneously. Typically, pilot sweeps with different phases or 
different frequency ranges are used to drive the different vibrators. The data 
are then correlated with each of the individual pilot sweeps to separate the 
data. Multiple sweeps are used to increase energy, and the cross correlations 
are added (stacked). Phase rotation of the sweeps also may be used to 
reduce hannonics. Separation of the data Is imperfect. Instead of clean 
vibrator records, the separated records may contain residual energy from 
other vibrators operating simultaneously. 

[0008] An alternative approach for separating vibrator signals and 
eliminating harmonics is taken by the IHigh Fidelity Vibratory Seismic i\^6thod 
(HFVS) disclosed in U.S. patents 5,719.821 and 5,721,710 to Sallas, et al. In 
the HFVS method the recorded seismic data are not correlated with a pilot 
signal, but instead are inverted using measured vibrator signatures from each 
sweep and each vibrator. Because the measured signatures include 
hanmonlcs, the inversion of the conresponding records recovers those 
hamnonics in the processed data, and thereby does not result In additional 
noise in the data. Because correlation is not used, correlation side lobes do 
not exist as a potential problem. Furthermore, Inversion with a measured 
vibrator signature can reduce effects fre>m variable vibrator coupling with the 
earth. However, in this method the vibrator motion for each data record is 
measured and used in the processing steps. The method includes use of a 
matrix inversion method to separate the signals from individual vibrators 
recording simultaneously. The matrix inversion requires that the number of 
sweeps M be greater than or equal to the number of vibrators N in order to 
solve a set of linear equations for the N vibrator signals. The ability to 
separate vibrator responses requires tiiat any two vibrators must differ in at 
least one of tiieir M sweeps. An advantageous way to accomplish this Is to 
phase-encode the M sweeps, typically witii one vibrator at a time sweeping 
with a phase shift relative to the other vibrators. The A/xiy^ vibrator signatures 
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are used to design a filter matrix that converts the M data records into N 
output records, one per each vibrator. Separation of vibrator records up to 60 
dB has been achieved with no visible degradation of the records from 
simultaneous recordings. 

5 [0009] The HFVS method is more fully described in association with 
Fig. 1, which depicts a typical land-based data acquisition system geometry, 
and Fig. 2. which depicts typical sweeps for four vibrators which may be used 
in that data acquisition system. Fig. 1 shows four vibrators 18. 20. 22. and 24. 
mounted on vehicles 34, 36, 38. and 40. The four different signatures 

10 transmitted into the ground during sweep i may be called sn, Si2, s,^, s,^. 
Each signature Is convolved with a different earth reflectivity sequence ^. ej. 
^ , 04 which includes reflections 26 from the interface 28 between earth layera 
having different Impedances (the product of the density of the medium and the 
velocity of propagation of acoustic waves In the medium). A trace <f/ recorded 

15 at a geophone 30 is a sum of the signature-filtered earth reflectivities for each 
vibrator. Formulating this mathematically, data trace d,(t) recorded for 
sweep i is: 



where ^^(0 = sweep i from vibratory. = earth reflectivity seen by vibrator j 
20 and 0 denotes the convolution operator. 

[0009] Persons sidiled In the art will understand the convolution 
operation and the convolution model upon which Equation (1) is based. Other 
readers may refer to standard treatises such as the Encyclopedic Dictionary 
of Exploration Geophysics, by R.E. Sheriff, 4* Ed. (2002), published by the 
25 Society of Exploration Geophysicists. (See the definitions of "convolution" 
and "convolutional model.") The noise term in Sheriff's definition of 
"convoiutionai model" has been neglected in Equation (1). This model is a 
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consequence of the concept that each reflected seismic wave causes its own 
effect at each geophone. independent of what other waves are affecting the 
geophone, and that the geophone response is simpiy the sum (linear 
superposition) of the effects of ail the waves. 

5 [0010] Thus, in this method vibrators radiate N sweeps into the 
earth, resulting In M recorded data traces. The HFVS method involves finding 
an operator, by solving a set of linear equations based on the known MxN 
vibrator signatures, that finds the set of iV' earth reflectivities that best predicts 
the recorded data. In the frequency domain, i.e., after Fourier transfonnation, 

10 the set of equations represented by Equation (1 ) are linear and can be written: 

I>iW) = t,S„(f)Ej(f) (2) 

7-1 

or, in matrix form for M sweeps and iV vibrators, 
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15 SE^D (4) 

If the number of sweeps is equal to the number of vibrators, this system of 
simultaneous equations can be solved for E : 



20 J = FD 

where 



(5) 
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F^(S)-\ (6) 

F is the filter or operator which when applied to the data separates it into 
individual vibrator records. 

[0011] For M ^ A/, Equation (4) may be solved by the method of least 
5 squares. For this mora general situation. Equation (4) can be written 

S'SE = S'D (7) 

where S' is the conjugate transpose of matrix S. Then, 

E = (s'sys'D, (8) 
and the filter F becomes 

10 F^(s*sys'. (9) 

[0012] The HFVS method can be used to record multiple source points 
simultaneously using a number of vibrators, but the use of more vibrators 
requires more individual sweeps each with its own listening time. It was 
previously not thought possible to eliminate the listening time, because the M 

15 sweep records must be separate unrelated measurements In order to solve 
the set of linear equations, i.e., otherwise the M equations (Involving N 
unknowns) would not be Independent. If the sweeps are cascaded without a 
listening time, then the reflection data from one segment would Interfere with 
data from the subsequent sweep segment, in addition, there would not be a 

20 one-to-one con^espondence between the data and the measured vibrator 
motions which represent the signatures put into the ground, so that hannonics 
would not be handled properly. The present invention solves these problems. 
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SUMMARY 

[0013] In one embodiment, the present Invention Is a method of 
operating N multiple seismic vibrators simultaneously at the same or different 
source points using continuous sweeps (i.e., realizing the benefits of multiple 
S sweeps without the disadvantage of a listen time at the end of every sweep), 
and separating the seismic response for each vibrator, which comprises the 
steps of (a) loading each vibrator with a unique continuous pilot sweep 
consisting of N segments; (b) activating the vibrators and using an anray of 
detectors to detect and record the combined seismic response signals from all 

10 vibrators: (c) recording the measured motion for each vibrator during the 
sweep; (d) parsing each motion record into N shorter records coinciding with 
the sweep segments, and then padding the end of each shorter record 
sufficiently to provide, in effect, a listening time for the shorter motion record; 
(e) forming an N x N matrix s whose element Sg{t) is the padded shorter 

15 vibrator motion record for the i* vibrator and /* sweep segment; (f) parsing the 
seismic data record from step (b) Into N shorter records, each shorter record 
coinciding In time with a padded shorter vibrator motion record; (g) forming a 
vector d of length N whose element d, is the i^' shorter data record from the 
preceding step; (h) solving for earth response Ejif) using the system of N 

20 linear equations and N unluiowns, SE = D, where Sf,(f) is the Fourier 
transfomi of Sy{t) and D,{/) is the Fourier transform of d,(t), where i = 1, 2. . 
. . N and J - 1, 2. ... N; and (h) inverse Fourier transforming the earth 
response Ej{f) back from the frequency (/) domain to the time (t) domain to 
yield e/t). 

25 [0014] In other embodiments, the vibrator motion records and the 
seismic data records may be parsed into M shorter records where M>N. In 
addition, the motions of the vibrators may be approximated by using 
theoretical wavefomis for the source signatures Instead of measured motion 
records. Harmonic noise is further reduced in preferred embodiments by 
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using sweep signals for the multiple vibrators that differ fronn each other only 
in the phase rotation of one or more sweep segments. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0015] The present Invention and Its advantages will be better 
5 understood by referring to the following detailed description and the attached 
drawings in which: 

[0016] Fig. 1 is a schematic diagram of a typical land-based acquisition 
geometry for four vibrators; 

[0017] Fig. 2 shows typical sweeps and listening time used to acquire 
10 i-IFVS data with four vibrators; 

[0018] Fig. 3 shows sweeps for four vibrators in one embodiment of the 
present invention; 

[0019] Fig. 4 shows sweeps for four vibrators in another emt)odiment of 
the present invention; 

15 [0020] Fig. 5 is a flow chart showing the main steps of one embodiment 
of the present invention; 

[0021] Fig. 6 shows the results of a model simulation of IHFVS 
recording; 

[0022] Fig. 7 shows the results of a model simulation of the present 
20 invention; and 

[0023] Fig. 8 shows the results of Fig. 7 for times longer than the 8-s 
sweep segment length. 

[0024] Changes and modifications in the specifically described 
embodiments can be carried out without departing from the scope of the 
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invention, wiiicii is intended to be limited oniy by the scope of the appended 
ciaims. 

DESCRIPTION OF THE INVENTION 

[0025] The present invention is a method for improving the efficiency of 
5 acquiring vibratory data with IHFVS techniques. With the HFVS method, data 
from a number of vibrators shaidng simuitaneousiy in seismic proximity to one 
another are separated by using (in one embodiment) a number of phase- 
encoded sweeps, where the number of sweeps is greater than or equal to the 
number of vibrators, resulting in a set of linear equations that can be solved 
10 simultaneously. The record length for each sweep includes an associated 
listen time containing reflections. The present invention eliminates the 
unproductive listening time for multiple sweeps but still provides the ability to 
separate the vibrator records and reduce contamination from harmonics. 
Production rates can be increased by as much as 30-80 %. 

15 [0026] In the present invention, multiple vibratory sources are used to 
record a land or marine seismic survey, and the signals are recorded by one 
or more detectors as shown in Fig. 1 for a land survey. Each vibrator Is 
excited by a different continuous sweep consisting of M segments, where M 
must be greater than or equal to the number of vibrators N, as iiiustrated in 

20 Fig. 3 for four vibrators and four segments, in the embodiment depicted in 
Fig. 3, each sweep segment is composed of a full sweep from conventional 
IHI=VS, eliminating the listening time between sweeps. Using such a sweep 
design, the vibrator begins and ends each sweep segment at rest, i-iowever, 
such a constraint is not necessary for the present Inventive method. 

25 Segmenting the sweep for each vibrator is a fundamental requirement of the 
present invention, but the sweep need not be designed around a pre-selected 
segment. Instead, segmenting may be an arbitrary, after-the-fact step of 
parsing a longer sweep into shorter parts, with the only requirements on the 
longer sweep being the same as for any vibrator sweep in conventional 
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vibroseis: the sweep should not exceed the inertial limitations of the vibrator, 
and the sweep should contain the full range of firequencles needed for target 
penetration and resolution. Two particular types of sweeps that may be 
familiar to vibrator users are upsweeps (steadily increasing frequency) and 
5 downsweeps (steadily decreasing frequency). Both work well in the present 
invention; e.g., the segments may be upsweeps or downsweeps. However, 
the present invention will work with segments that employ any physically 
realizable sweep type Including linear, nonlinear and pseudo-random. In 
some preferred embodiments, each sweep segment is a linear up or down 
10 sweep that encompasses the full range of frequencies required for imaging 
the reflectors with the desired resolution as would be familiar to those who 
design seismic acquisition surveys. In other embodiments, the sweep 
segments could be composed of a combination of upward and downward 
progressions in frequency, nonlinear sweeps, or pseudo-random sequences. 

15 [0027] In prefen'ed embodiments of the present invention, the duration 

of each segment is greater than the travel time for the target reflectors or 
greater than the listening time for conventional recording. Just as with the 
correlation process, inversion of vibrator data compresses the energy from a 
sweep or sweep segment into a pulse. If a second sweep segment begins 

20 sooner than the travel time to and back from the deepest reflector of interest 
(the "listening time"), then shallow reflections from the second segment can 
interfere with the deep reflections from the first segment. A dead time of any 
length can be inserted between segments, but preferably the dead time is 
zero because this is the value of the invention. As used herein in connection 

25 with the present invention, a "continuous sweep" means a sweep in which any 
dead time between sweep segments is preferably zero, but never more than 
the listen time used in conventional HFVS. 

[0028] The sweep for each vibrator must be unique, i.e., no two 
vibrators can have exactly the same sweep. This is necessary in order that 
30 the data for each vibrator can be separated from the other vibrators. One way 
of accomplishing this is by applying a phase rotation to one or more of the 
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segments of the sweep. For example, the first vibrator can have the first 
segment of the sweep at 90 degrees from all the other segments. The 
second vibrator can have the second segment at 90** from the other 
segments. The third vibrator has the third segment at 90* from the other 

5 segments. This pattern can continue for N vibrators and N segments as 
shown in Fig. 3. Alternatively the phase angles discussed in the HFVS patents 
can be used for the different segments. Using the cascaded sequence 
described by Anderson in which subsequent segments are phase rotated by 
0. 90, 180 and 270 degrees combined with the 90" phase rotation above, as 

10 shown in Fig. 4, has particular advantages in reducing hamnonlcs. Other 
ways of generating a unique sweep by using different frequency ranges, 
sweep rates, or by using different random sweeps can also be used in the 
present invention. 

[0029] Fig. 5 is a flow chart showing the main steps of one embodiment 

IS of the present invention. In step 101, the sweeps, as shown for example in 
Figs. 3 or 4, are loaded into the vibrator controllers for the conesponding 
vibrators. Each vibrator receives a unique sweep, conveniently (but not 
necessarily) accomplished by the phase rotation technique discussed 
previously. Unlike the HFVS method, the sweep is a single, long, continuous 

20 sweep, exemplified by the sweeps shown in Fig. 3 and Fig. 4 which are 
composed of four segments of 8 j each, followed by an 8 j listen time. This is 
to be contrasted with the typical HFVS sweep of Fig. 2 in which there are four 
8 s sweeps, but each is followed by an 8 j listen time. Thus, in this example, 
to get the same amount of source energy into the ground takes 64 s with 

25 HFVS as compared to 40 s with the present inventive method. Each of the 
single continuous sweeps used in step 101 for the i\r vibrators must be divided 
into at least N segments, each segment in the sweep of any one vibrator 
being of the same length (time duration) as the corresponding segments In 
the sweeps of all the other vibrators. For example, the second segment 

30 should be the same length for all vibrators, as should the fourth (or any other) 
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segment, but the length of the second segments can be different than the 
length of the fcMjrth. 

[0030] In step 102 (could be perfonned before step 101), the vibrators 
are located at preselected locations. All the vibrators are then simultaneously 

S exdted by their conresponding pilot sweep (step 103), and a single long 
record is recorded from one or more detectors in step 104. The length of the 
data record will be the length of the pilot sweep plus one listening time. In 
addition, the measured motions of each vibrator, typically signals from 
accelerometers mounted on the baseplate and on the reaction mass of each 

10 vibrator are recorded. Besides its use in the s and S matrices in Equations (1)- 
(5), this ground force signal, which may be computed as the mass-weighted 
sum of the baseplate and reaction mass accelerometer signals, is typically 
used in a feedbacl^ loop to control the excitation of the vibrator. Any other 
signal such as the pilot signal Itself that could be considered representative of 

IS the source signature may be used for the purposes of the present invention. 

[0031] in step 105, the measured motion records for each vibrator are 
parsed Into M traces (shorter records) composed of iheM^N individual time 
segments into which the pilot sweeps were divided in step 101. The shorter 
records are then lengthened by adding zeros to the end, called padding the 
traces. The end padding should extend the duration sufficiently to, in effect 
(not in reality), provide a listening time (two-way seismic wave travel time to 
the deepest reflector of interest) for the segment. More padding beyond that 
desired amount will increase computation time without added benefit. Zeros 
or padding can also be applied to the beginning of each trace if desired. 
Because the segment length is preferably chosen to be longer than the 
desired listening time, it may be convenient to standardize the total padding 
duration to be the same as the original segment duration so that the total 
trace length is twice the segment length. The M traces become part of the 
sweep matrix s from Equation (1). For example, s^^ is equal to the measured 
motion for the first vibrator for the duration of the first segment plus the zero 
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padding at the end of the segment. The element s^2 Is equal to the measured 
motion for the first vibrator for the duration of the second segment plus the 
zero padding at the end of the segment, etc. 

[0032] In step 106, a copy of the geophone data record Is parsed to 
5 make N shorter recorels, each of duration equal to a segment duration plus the 
duration of the padded length used for the measured motions. In one 
embodiment of the present invention, the first record would con^espond to the 
first segment plus the data before and after the segment conresponding to the 
padded time. The second record would consist of the second segment plus 
10 the padded time, etc. The N records mal<e up the vector d that appears in 
Equation (1) above. The end padding in step 105 should be of sufficient 
duration to capture the seismic response due to the end of that source motion 
segment. 

[0033] If desired, more sweeps can be performed to build up energy, 

15 adding more rows to the sweep matrix s and the data vector d . Because the 
system of simultaneous equations will not be linear in the time domain, the 
Fourier transform is computed yielding the matrix S and vector D, and a 
separation and inversion filter F is derived In step 107 by inverting the matrix S 
using Equation (9). Equation (6) may be used if Mis chosen to be equal to N. 

20 For Af > N, the system of equations Is over-detennined, and a best-fit solution 
is obtained using a criterion such as least squares. This approach may be 
useful even where M is chosen to be equal to N because one of the vibrators 
may be temporarily unavailable and recording could proceed with fewer 
vibrators. The method can proceed even if the number of vibrators drops 

25 temporarily to one during acquisition. A person skilled in the art will 
understand that although matrix fomnalism is used in the preceding 
description, any method, numerical or analytical, of solving M simultaneous 
linear equations in N unknowns may be used in the present invention. Next, 
in step 108, the filter is applied to the data vector D, and the inverse Fourier 

30 transform is computed resulting in N separated records Ej0 where / is 
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frequency. With the above procedure, harmonics and data from subsequent 
and previous segments will appear at times greater than the segment length, 
and they will not Interfere with the target reflections. Persons skilled in the art 
will understand that the filter F could be inverse-transformed to the time 
5 domain and then applied to the time domain data. Similarly, the order of 
steps 106 and 107 may be Interchanged. Such aKemative procedures are 
insubstantial changes and hence equivalents to the procedure described 
above, and therefore are part of the present invention. 

Examples 

10 [0034] Fig. 6 and Fig. 7 compare the firet 2.5 seconds of model data 
after the process of separation and inversion using the conventional HFVS 
method (Fig. 6) and the method of the present invention (Fig. 7). The model 
data are generated using 51 receivers. 400 feet apart. Four sources are 
located at distances of 5000. 8333, 11666. and 15000 feet from the first 

15 receiver. The data for each source location were convolved with actual 
vibrator signatures from field measurements and combined to simulate 
simultaneous acquisition of the four sources. The vibrator sweep for HFVS 
was an 8-j linear sweep from 8 to 128 Hz. The vibrator sweep for the present 
Invention is a Z2-s sweep composed of four 8-j segments. The horizontal axis 

20 in each figure displays receiver location, by source. The results after 
separation and inversion of the model data show little difference between the 
two methods, and both perfectly separate the reflection data for the 4 sources. 

[0035] Fig. 8 shows the separated and inverted results generated by 
the present invention for a longer time period of the continuous sweep than is 

25 shown in Fig. 7. The sweep segment Is a downsweep. As can be seen, the 
separated records are clean for the top 8 seconds, which corresponds to the 
segment length as illustrated in Fig. 3. (This Illustrates the reason why 
segment length in the present invention is preferably chosen to be at least as 
long as the seismic wave travel time down to and back up from the deepest 

30 reflector of interest.) Interference from the subsequent sweep occurs after 8 s 
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and interference from the previous sweep occurs after ^5 s. These are 
separated into the individual source locations. Noise from harmonics, which 
are not perfectly matched with the proper vibrator signatures, appears after 
the primary interference from the subsequent sweep at around 10-14 
5 seconds. 

[0036] The foregoing description is directed to particular embodiments 
of the present invention for the purpose of illustrating it. It will be apparent, 
however, to one skilled in the art that many modifications and variations to the 
embodiments described herein are possible. For example, In step 107 of Fig. 

10 5, the system of equations Is Fourier transfonned to the frequency domain. 
Any other transform that produces a system of linear equations in the 
transform domain will work in the present inventive method, and the 
appended claims are to be understood to include any such transforms. All 
such modifications and variations are intended to be within the scope of the 

15 present invention, as defined in the appended claims. 



